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O ■ Abstract 
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■ We present an explicit exact solution of Einstein's equations for an inhomo- 

geneous dust universe with cylindrical symmetry. The spacetime is extremely 
simple but nonetheless it has new surprising features. The universe is "closed" 
in the sense that the dust expands from a big-bang singularity but recollapses 
to a big-crunch singularity. In fact, both singularities are connected so that the 
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(3JT), whole spacetime is "enclosed" within a single singularity of general character. 

The big-bang is not simultaneous for the dust, and in fact the age of the universe 
as measured by the dust particles depends on the spatial position, an effect due 
to the inhomogeneity, and their total lifetime has no non-zero lower limit. Part 
of the big-crunch singularity is naked. The metric depends on a parameter and 
contains flat spacetime as a non-singular particular case. For appropriate values 
of the parameter the spacetime is a small perturbation of Minkowski spacetime. 
This seems to indicate that flat spacetime may be unstable against some global 
non-vacuum perturbations. 

Exact solutions of Einstein's equations for an inhomogeneous dust with an Abelian 
(j2 on S*2 group of motions seem to be difficult to find, and very few of them are known, 
see [p] and references therein. Among them, there are some Petrov type-D known cases 
belonging to the general class found by Szekeres 0, some examples appear in J3, ffl. 
However, it seems that no exact solution for an inhomogeneous dust with cylindrical 
symmetry and references therein) has been singled out. In this short letter we prove 
that, actually, there is a very simple cylindrically symmetric dust universe included in 
the Szekeres family, which surprisingly has not been studied or considered hitherto. 
Some plausible reasons for this overlooking will also be analized: apparently, some 
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minor errors in led to a misunderstanding which has kept the solution "behind the 
screen" . 

As a matter of fact, the solution was recently rederived as the only dust solution 
defined by the following separation Ansatz in half-null coordinates ||: 

ds 2 = -2dudv + [f{u) + g{v)] 2 dy 2 + [h{u) + k(v)f dz 2 , 

where none of the derivatives of the given functions vanish, and the coordinates are 
taken to be comoving, that is, the dust velocity vector field reads u oc d/du + d/dv. 
Changing to typical Lorentzian coordinates the line-element takes the strikingly simple 
form 

ds 2 = _ dt 2 + df) 2 + p 2 d(f 2 t±f^ (1) 

where a is an arbitrary constant. The cylindrical symmetry is explicit, where p > 
is the cylindrical radius (the axis is at p = 0) and <p is the corresponding angular 
coordinate running from to 2tt. It can be checked that there are no further isometries. 
The coordinate z is taken to run from — oo to oo, while the coordinate t will be restricted 
by the singularity of the spacetime, as we will see presently. 

The dust is comoving in the coordinates of (|I|), its energy density given by 

e = 4 (a 2 - t 2 - P 2 ) 1 , 

while the non- vanishing components of the kinematical quantities for u — d/dt com- 
puted in the orthonormal co-basis given by 6 a oc dx a are 

6 = -2t (a 2 -t 2 -p 2 )' 1 , 
6 J 

011 = 022 = —; CT 33 = 2-. 

Of course, the acceleration and vorticity vanish. The scalars of the Weyl tensor, 
computed in the null tetrad (see 0) fe = 2~ 1 / 2 (0° - 1 ), I = 2~ 1 / 2 (0° + 1 ), m = 
2- 1 / 2 (0 2 + 20 3 ), read 

^0 = ^4 = 3^ 2 = -~Q, i = ^ 3 = 0, 

and thus these solutions are of Petrov type D. 

The expressions of the energy density and of the Weyl tensor show that a curvature 
singularity appears at t 2 +p 2 — a 2 — > 0, which can be viewed as a connected hypersurface. 
The meaningful spacetime region (defined by g > 0) is clearly given by 

t 2 + p 2 < a 2 , 

and therefore the singularity 'wraps' the entire manifold. This singularity has a general 
character, as it is spacelike in the regions given by t E [—a, —a/\^2) U (a/v2, a] (so 
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Figure 1: This diagram corresponds to the (t,p) plane of the spacetime with line- 
element given by ([!]). As usual, null lines are at 45°. The whole spacetime is the 
product of this plane with the (p, z)-surfaces, which correspond to the orbits of the 
2-dimensional isometry group defining the cylindrical symmetry The thick curve rep- 
resents the spacetime singularity. The dust particles move along the geodesies repre- 
sented by arrowed vertical lines in the physical region where g > 0, whose "size" is 
measured by the radius a of the semi-circle. In the limit with a — > oo the singularity 
disappears and the solution becomes flat Minkowski spacetime. 

that p G [0, a/y/2)), null in t — ±a/y/2 (where p = a/y/2), and timelike at the region 
* G (-a/y/2, a/V2) (thus p G [a I v^2, a]). This is shown in Figure [I]. 

All the causal geodesies start and end at the singularity. As a particular case, the 
timelike geodesies defined by the dust flow begin at the part of the singularity where 
t < with 9 — > oo, go on expanding up to t = 0, where the dust starts to contract and 
eventually die in the future part of the singularity (t > 0) where 9 — > — oo. Actually, 
all endless causal curves passing through any given point of the spacetime reach the 
singularity both in their past and future for a finite value of their generalized affine 
parameter. Thus, the singularity is past and future universal and can be termed as a 
big-bang/big-crunch singularity with a general character (see || for definitions). 

Some interesting features can be deduced from the diagram of Figure [I]. To start 
with, the big-bang and the big-crunch are not simultaneous, something which happens 
in some other well-known cases, see e.g. |§. The total proper lifetime of any dust 
particle is always finite and depends on its position (on p), so that it is smaller for 
bigger p. This total proper time has no lower bound, so that one can find dust particles 
with as short a lifetime as desired. On the other hand, the maximun lifetime for a dust 
particle is reached at the axis, and is given by t max = 2a. This increases with a. Notice 
that for very big a the metric approaches flat spacetime, and in the limit a — > oo the 
spacetime is exactly Minkowskian (and the singularity disappears!). 
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Actually, for small values of 1/a the metric can be considered as a small pertur- 
bation of the flat Minkowski metric. As a is a free parameter, the smallness of the 
perturbation can be chosen at will. Notice also that the spacetime is obviously globally 
hyperbolic and, furthermore, the hypersurface t = is maximal in the sense that its 
second fundamental form vanishes. Thus, this spacetime seem to indicate that the clas- 
sical results on the global non-linear stability of Minkowski spacetime [|K], may not 
be generalizable to the case of non-vacuum perturbations (or to non-asymptotically-flat 
perturbations). It is interesting to see how this metric shows that the appearance of 
an arbitrarily small quantity of matter everywhere in flat spacetime destroys its global 
structure and encloses the matter within a strong curvature singularity. 

Another curious property is that there are causal curves (and geodesies) which 
never cross the t = hypersurface, starting and ending at the part of the singularity 
with t < 0. Thus, a part of the singularity with t < is actually the future singularity 
for some physical particles (not for the dust, though). In fact, some dust particles 
(and other particles as well) can communicate with each other at their respective birth 
times. In other words, some dust particles can send causal signals (even from their 
bang starting time) which reach other dust particles at their corresponding birth times. 
Thus, for instance, a dust particle could tell another dust particle (with a big p) about 
its fate, so that the latter may decide to leave the dust motion and travel to smaller 
values of p in order to live longer. Similarly, a dust particle dying in the big-crunch 
can send a signal at that very moment showing its fate to other dust particles which 
could see how their end looks like. In other words, the big-crunch singularity is partly 
"naked" and visible for the dust. In general, and as the metric is time-symmetric, all 
statements for the future have their corresponding past counterpart. 

The line-element ([I]) is a dust solution with Petrov type D, so that it belongs to 
the family of the so-called Szekeres spacetimes 0, || [12| , which were later generalized 
to a bigger family of Petrov D spacetimes [T^, [TJ], |S| known as the Szafron-Szekeres 



solutions. The Szekeres solutions do not contain any isometry in general |Ld[, but 
particular families of this kind of inhomogeneous solutions with a G2 group of isometries 
were given in || and 0] . More precisely, the metric (|lj) is in fact contained in the case 
PII of the parabolic models in the classification of the Szekeres solutions presented 
by Bonnor and Tomimura in || (page 86) where the evolution of these models was 
studied. 

Nevertheless, this solution was not under the scope of the aforementioned references 
mainly due to the fact that interest was focused on the study of the solutions containing 
only an initial spacelike singularity. Indeed, in 0] it was already established that the 
Szekeres solutions contain two possible singularities, corresponding to the vanishing 
of two functions appearing in the metric: the so-called Friedmann function R(t) and 
another arbitrary function Q(t, x, y, z) (in what follows we are using the notation of [[|). 
The singularity coming from R = occurs on a spacelike hypersurface t = 0, whereas 
that from Q = admit a great variety of behaviours, and thus in ||] the assumption 
that t = were reached before Q = was tacitly demanded. The arbitrary functions 
ought then to be chosen such that Q > in the spacetime in order to have g > 0. 
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However, in the PII class one has R = 1 so that the only singularity corresponds to 
Q = 0. As we have seen, the line-element (P explicitly reaches Q = 0, which now is 
simply Q = 1 — (t 2 + p 2 )/a 2 . The rest of the functions as they appear in || are given 
as follows: fi = v = a = 0, u = 1, (3 = const. (= — 1/a 2 ). In fact, it was argued in 
H that the PII class would always contain a region of the spacetime with a negative 
energy density, but this is not the case as we have explicitly proved for the solution 
under consideration: for the solution ([I]) the condition Q > is enough to define a 
spacetime with g > everywhere. 

Regarding [IB], the class PII was again dismissed because it has no non- vacuum 



Friedmann-Lemaitre-Robertson- Walker (FLRW) specialization, so that its singularity 
structure, as well as the study of its evolution, cannot be performed with a reformulated 
line-element in terms of increasing and decreasing density perturbation modes of FLRW 
dust models, which was the technique used. 
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